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1. Introductory Remarks and Calculation Setup 

Flavour Changing Neutral Currents (FCNC) and CP violation may furnish useful informa- 
tion on the impact of models defined beyond the standard model (BSM). In various BSM models 
(like for example the supersymmetric ones) there appears the possibility for AS = 2 processes at 
one loop, even mediated by the strong interactions. These effects are thus potentially large. The 
computation of the relevant matrix elements of the effective Hamiltonian in combination with the 
experimental value of Ek offers the chance to constrain the values of the model parameters (like 
for instance the off-diagonal terms of the squark mass matrix in supersymmetric models [jl|]) which 
enter explicitly in the Wilson coefficients. 

In all the BSM models the effective Hamiltonian relevant for the AS = 2 processes takes the 
general form 

=^ef = (1-1) 

where the operators i^,- are defined by 

^1 = Wy^{l-Ys)d"][^y^{l-y,)d% 

^2 = [^'{\-Ys)d"][^{\-Y5)d% i^, = [f{\-Y,)d%i{\-Ys)d"l 

= r{\-Y5)d"][s'{l + Y5)d% ff5 = [s"{l-Y5)d'][s'{l + Y5)d"] (1.2) 

We note that in the SM case only the operator contributes. The parity-even parts of the operators 

^1 = [s"Y,{^ + Y5)d"][s'Y,i^ + Y5)d% 

^2 = [s"{l + Y5)d''][s'{l + Y5)d% ^s = [s"{l + Y5)d''][s'{l + Y5)d'] (1-3) 

coincide with those of the operators i^,. Therefore, due the parity conservation of the strong in- 
teractions only the parity-even contributions of the operators need to be calculated. Defining a 
basis of the parity even operators as follows 

O^^ = {sY^d) {sYnd) , = {sYt, Ysd) (^"Tm Ysd) , 

(9^^ = {sYsd) (sYsd) , O^^ = {sd) (sd) , 

O^^ = {so^yd){so^yd) (1.4) 
through a Fierz transformation we obtain 

^4 = (O^^-O^^), ^5 = _l(o^^_o^) (1.5) 

Up to now, lattice calculations have been presented in the quenched approximation ([^, [^], 
[Q]) with the exception of a preliminary study of the bare matrix elements using unquenched sim- 
ulations with 2-1-1 dynamical quarks |^. 

Our lattice computations have been performed at three values of the lattice spacing using the 
Nf = 2 dynamical quark configurations produced by the ETM collaboration [|6|]. ETMC dynamical 
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configurations have been produced with the tree-level Symmanzik improved action in the gauge 
sector while the dynamical quarks have been regularized by employing the twisted mass (tm) for- 
malism 1^. It has been demonstrated that with the condition at maximal twist this formalism 
provides automatic 0(a)-improved physical quantities [||]. 

In the so called physical basis the fermion lattice action concerning the sea sector is written 

5j^a'° = «'Iv^W{^l7A,(VM + v;)-/75T^[Mcr-^IV*V^]+Ms.«WW (1-6) 

where the Wilson's r parameter has been set to unity, ^lf{x) is the quark flavour doublet, and V* 
are nearest-neighbour forward and backward lattice covariant derivatives, pLsea is the (twisted) sea 
quark mass and Mcr the critical mass. It has been shown that the use of the tm regularization can 
simplify the renormalization pattern properties of the four-fermion operators (e.g. Bk) [0> H Hi- 
Moreover, both 0{a) improvement and continuum-like operator renormalization pattern can be 



achieved introducing a valence quark action of the Osterwalder-Seiler type [ 1 1 1 by allowing for a 



replica of the down (<i, d') and strange {s, s') flavours [12]. The valence quark action assumes the 
form 



Svai = I I ^/ W { ^ I 7m (Vm + - ijsrf [M„ - ^ £ V* V^] + ^/j ^/ W (1-7) 

X f=d,d',s,s' H 

with — = rd = r^' = r^i = 1. Note that the field qf represents just one individual flavour. The four 



fermion operators of Eq. (1.4) can be written in general form as Oif = 2{[^ir^2]fer^4] + {l2 ^ 
^4)} with qi and ^3 identified with the strange quark (by setting /i^ = /ly/ = jJ-st range) and q2 and ^4 
identified with the down quark (by setting /i^ = /i^/ = /i^); the interpolating fields for the external 
(anti)Kaon states are made up of a tm-quark pair {dy^s, with — = rj) and a OS-quark pair (d'y^s', 
with r^i = r,/). This mixed action setup with maximally twisted Wilson-like quarks has been studied 



in detail in Ref. [ ]12p and it has been demonstrated that it allows for an easy matching of sea and 
valence quark masses and leads to unitarity violations that vanish as a^ as the continuum limit is 
approached. Moreover in the present computation the quark mass matching is incomplete because 
we are neglecting the sea strange quark (i.e. we work in a partially quenched set-up). A first test 
that the proposed method leads to 0{a) improved results was already performed in the calculation 
of Bk with fully quenched quarks [|^]. In a recent publication [0] our collaboration, using the OS- 
tm mixed action set-up, has presented an 0(a) -improved computation of Bk with Nf = 2 dynamical 
quarks. Using non-perturbative operator renormalisation and three values for the lattice spacing, 
the RGI value of Bk in the continuum hmit is = 0.729 ± 0.030. 

In Table [l| we give the simulation details and the values of the sea and the valence quark 
masses at each value of the gauge coupling for the calculation presented in this work. The smallest 
sea quark mass corresponds to a pion of about 280 MeV for the case of j8 = 3.90. For j8 = 4.05 
the lightest pion weighs 300 MeV while for jS = 3.80 the lowest pion mass is around 400 MeV. 
The largest sea quark mass for the three values of the lattice spacing is about half the strange quark 
mass. For the inversions in the valence sector we have made use of the stochastic method (one-end 
trick of Ref. [p3|]) in order to increase the statistical information. Propagator sources have been 



located at randomly chosen timeslices. For more details on the dynamical configurations and the 



stochastic method application see Ref. [16] 
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0.065 fm) 
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Table 1: Simulation details 
2. B-Parameters and Four-Fermion Matrix Elements 



As it has been shown in [|12|], the discrete symmetries guarantee that in the OS-tm mixed 
action set-up the renormalisation of the four-fermion operators is continuum-like in the sense that 
the mixing between operators of different naive chirality is of order 0{a^) or higher. An equivalent 
view of the same property can be offered by the fact that in the (unphysical) tm-basis the parity- 
even part of each of the four fermion operators is mapped over its parity-odd counterpart. Then 
due to the CPS symmetries the parity odd operators have the same block-diagonal renormalisation 



matrix pattern both in the continuum and at finite value of the lattice spacing ([|17[], []1 8[]) 
The B -parameters for the operators dOb are defined as 



(^|^i(/l)|^:") = B,{^l)-mlfl^BK{^l)-myi 

where Q = {—5/3, 1/3,2,2/3}, / = 2,...,5. The matrix element of the operator vanishes in 
the chiral limit while the matrix element of the operators / = 2, . . . , 5 get a non-zero value in the 
chiral limit. From the above equations it can be seen that the calculation of the B, parameters for 
/ = 2,...,5 involves the calculation of the quark mass at the same scale /i. In order to avoid 
any extra systematic uncertainties in the computation of the matrix elements due to the quark 
mass evaluation, it has been proposed the calculation of appropriate ratios of the four-fermion 
matrix elements ([@], Here, besides the calculation of the B parameters, we also consider the 
following ratios 

R-(it\ « {K'\m)\K') ^ 

The computation of the renormalisation constants (RCs) relevant for both the four-fermion and 
two-fermion operators' has been performed in a non-perturbative way using the RI-MOM scheme 
following the strategies detailed in Refs. ^ and ||^]. In Fig. [I] we show, for jS = 3.90, that all the 

'in our OS-tm mixed action set-up we need to use the RCs for the scalar and the pseudoscalar density operators in 
the calculation of B,- for i = 2, ... ,5. For Si, instead, (i.e. Bk) the normalisation constants for the axial and the vector 
current are needed. 
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Figure 1: The off-diagonal off-block elements of the four-fermion RC-matrix operator, A/j , (for j3 ~ 3.90) 
which take values compatible with zero. 



off-diagonal off-block elements of the four-fermion RC-matrix operator, A,y , which are expected to 
vanish with the tm-OS mixed action setup, take values compatible with zero. In Fig. ^ we show the 
combined fits of / = 2, . . . ,5 with respect to the renormalised light quark mass, ro/if , in 

the MS scheme. 

In Table we present our preliminary results in the continuum hmit and in the MS scheme for 



the B -parameters and the ratios 



<0i> 



/ = 2, . . . ,5 calculated at the physical point {}JLd,i^s)- The 



ratios, ^^j^ , / = 2, . . . ,5, have been calculated either directly (through /?,) or using the B; estimates 
and the values of the u/d and strange quark mass [23]. The results are compatible within one or 
two standard deviations. We have tried fit functions using either a second or first order polynomial 
with respect to the light quark mass to which a term proportional to has been added; we do not 
notice a significant difference in the final continuum limit values. We should note that the use of a 
fit function containing a NLO logarithmic term leads to rather similar results with those obtained 
with a second order polynomial fit function. 
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Figure 2: Combined fits of ; = 2, . . . ,5 with respect to the renormaHsed Hght quark mass, rQ^f, 

in the MS scheme. The fit functions, shown here, are second order polynomial functions of the renormalised 
light quark mass with the addition of a term proportional to a^. 



fit function / mi (using/?,) (using fiQ Bj 



quadratic 2 -15.5(1.0) -17.1(1.5) 0.56(0.04) 

3 7.1(0.5) 8.7(0.9) 1.43(0.13) 

4 24.6(1.4) 27.9(2.5) 0.76(0.06) 

5 6.9(0.7) 7.6(1.2) 0.63(0.09) 



linear 2 -15.0(0.6) -17.1(1.2) 0.56(0.02) 

3 7.0(0.3) 8.8(0.6) 1.44(0.08) 

4 24.2(0.9) 27.8(1.8) 0.76(0.04) 

5 6.6(0.5) 7.5(0.7) 0.62(0.06) 

Table 2: Preliminary results in the continuum limit for the B-parameters and the ratios , / = 2, . . . ,5. 
calculated at the physical point (/i^/, /ii, a — 0). All results are given in the MS scheme. The ratios of the 
operators' matrix elements are estimated using either the direct (3rd column) or and the indirect method i.e. 
through the B, calculation (4th column) . 
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